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'"^ , Abstract. We study a large deviation functional of density fluctuation by analyzing 

D ' stochastic non-linear diffusion equations driven by the difference between the densities 

fixed at the boundaries. By using a fundamental equality that yields the fluctuation 
theorem, we first relate the large deviation functional with a minimization problem. 
We then develop a perturbation method for solving the problem. In particular, by 
. performing an expansion with respect to the average current, we derive the lowest 

order expression for the deviation from the local equilibrium part. This expression 
implies that the deviation is written as the space-time integration of the excess entropy 
production rate during the most probable process of generating the fluctuation that 
corresponds to the argument of the large deviation functional. 

O 

PACS numbers: 05.10.Gg, 05.70.Ln, 05.10.-a 

> 

CO, 1. Introduction 

o 

The nature of macroscopic fluctuations in equilibrium systems is determined by 
\Q ' thermodynamics. Mathematically, this implies that the large deviation functional for 
the fluctuations is given by a thermodynamic function. This relation is called Einstein's 
formula, and it is an important consequence of equilibrium statistical mechanics. 

Apart from thermodynamic systems, there are many systems in which macroscopic 
^ ■ fluctuations are studied. In general, it appears sound to consider that properties of 
macroscopic fluctuations for non-equilibrium systems depend on the systems under 
consideration. However, from an optimistic viewpoint, it may be expected that 
there exists a class of non-equilibrium systems in which the large deviation functional 
has a physical correspondence with macroscopic deterministic laws such as those in 
thermodynamics and hydrodynamics. 

Recently, researches speculating the existence of such a system have been presented. 
In one approach, an extended form of thermodynamics is explored before considering 
fluctuations [lj. If it is operationally and consistently constructed, this extended 
function might provide a large deviation functional. Although there are some successful 
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examples [21 E] , it turned out a severe restriction is required to extend a thermodynamic 
framework in macroscopic systems [3]. 

Another approach utilized a new variation principle called the additivity principle 
[H El E]. Indeed, with the aid of this principle, non-trivial forms of large deviation 
functionals were exactly derived for some models. Although the discovery of this 
principle stimulates us to expect it to possess generality, there is no evidence that it is 
useful for a wide class of non-equilibrium systems. 

Considering the above-mentioned difficulties in constructing a general framework, 
we have decided that concrete expressions of large deviation functionals should be 
investigated for several models. With this motivation, we attempt to develop a 
perturbation method for calculating large deviation functionals. In particular, we 
analyze a fluctuating hydrodynamic model that describes the density fluctuations in 
systems that are attached to particle reservoirs at the boundaries (see Sec. 12). Even 
for this simple model, a perturbation method has never been proposed, except for the 
approximation method by which large deviation functionals are treated as a quadratic 
form [7J. This approximate expression can be obtained by the analysis of a linearized 
equation, but it is not appropriate for our purpose because we wish to know the physical 
correspondence of large deviation functionals. 

In this paper, we expand a large deviation functional with respect to the average 
current, maintaining the nature of large fluctuations. Here, it should be noted that 
the large deviation functional can be exactly derived for equilibrium cases. This exact 
derivation is due to a detailed balance property. Thus, in order to develop a perturbation 
method with respect to the average current, we need to find a useful relation that 
corresponds to an extended version of the detailed balance property. Fortunately, a 
fundamental relation that recovers the detailed balance condition in equilibrium cases 
is known to be related to the fluctuation theorem [H EJ [101 [HI [12] . Using this relation, 
in Sec. EJ we formulate a minimization problem to which a perturbation method can be 
applied for determining the large deviation functional. In Sec. HI we calculate the large 
deviation functional of the lowest order with respect to the average current. Section E] 
is devoted to concluding remarks. 

2. Preliminaries 



We consider a one-dimensional system in contact with particle reservoirs at the 
boundaries x = and x = L. Let p(x, t) be the coarse-grained density field. We 
assume that p(x, t) obeys the stochastic partial differential equation 



2.1. Model 




(1) 




(2) 
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where D in (pQ) and a in represent the diffusion constant and the conductivity in the 
system under consideration, respectively. T is the temperature, and the Boltzmann 
constant is set to unity. It is straightforward to extend our analysis for d(> 2)- 
dimensional systems, but for simplicity we focus on the one-dimensional case. 

Since we are interested in the large deviation functional of the density field, we 
set e = 1/L and introduce a new space-time coordinate (x',t') = (ex,e 2 t). With this 
coordinate, we define a new field p'(x',t') = p(x,t) and noise £ r (x r ,t') = e _1 ^(a;, t). 
Then, p'(x',t') satisfies the same equation as (CD) with prime symbols. Omitting the 
prime symbols, we obtain ([I]) with 



We impose the boundary conditions p(0, t) = po and p(l, t) = pi, where we assume that 
Po < Pi without loss of generality. 

2.2. Notation 

Since a few types of functional dependence appear in this paper, we use the following 
notations in order to avoid confusion. Formally, p represents the density field as a 
function of (x,t), where < x < 1 and < t < oo. However, as in the usual 
practice in physics, we write p(x, t) as p if the argument can be easily guessed. On 
the contrary, when a quantity g depends on the function p, we explicitly write the 
functional dependence as g(pQ). Thus, for example, -D(p) in ([1]) does not represent 
the functional dependence of p, but it should be read as D(p(x,t)). Furthermore, the 
density field as a function of x is denoted by p. In a similar manner, h(p) represents 
h(p(x)) and g(p()) represents the functional dependence of p. We also use the expression 
p( , 0) = pQ as an identity for functions of x. That is, this means that p(x, 0) = p(x) 
for any value of x. 

2.3. Question 

We study statistical properties in the steady state of this model. The average current 
J and average density profile p(x) are obtained from the relation 



with the boundary conditions p(0) = p and p(l) = p 1 . However, it has been known 
that fluctuations around the average value exhibit non-trivial behaviors such as non-local 
correlations [?]. Here, let P s (pQ) be the stationary distribution of density fluctuations. 
In the limit e — > (L —>■ oo in the original problem), it is characterized by the leading 
order expression 



(£(x, my, s)) = 2ea(p(x))TS(x - y)5(t - s). 



(3) 



D(p)d x p = - J 



(4) 



P s (p())^exp --I(pQ)) 



(5) 



where I(p()) is called the large deviation functional. 
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Despite the simplicity of the model, it is quite difficult to calculate the large 
deviation functional for general non-equilibrium cases. To date, only for the special 
case where D = 1 and a(p) = p(l — p), it was derived as a variational form by using the 
Hamilton- Jacobi formulation for the model described by ([1]) with ([3]) [13]. It should be 
noted that this large deviation functional is identical to the exact solution for a simple 
exclusion process in contact with particle reservoirs at the boundaries [I]. However, we 
find that the special technique for solving the Hamilton- Jacobi equation cannot be used 
for other functional forms of D(p) and cx(p). 

Here, let us consider the problem from a physical viewpoint. We first review the 
result for the equilibrium case po = Pi- According to Einstein's formula, the large 
deviation functional is written as 

I(P0) = / WO) - F(pQ)] (6) 

with 

F(pQ)= I dx[f(p(x))-w(x)), (7) 
Jo 

where (3 is the inverse temperature; f(p), the free energy density; p, the chemical 
potential defined by p — f'{po)- Hereafter, the prime represents the derivative with 
respect to the density. Furthermore, with regard to the intensity of density fluctuations 

X{— /q dx(p(x,t) — po) y /e), we know that the equality /" = Tx _1 and the 

Einstein relation D\ = o~T hold (see the review in Ref. [H]). These lead to 

In this manner, we have determined the large deviation functional without any 
calculation. On the contrary, without any physical consideration, we can derive the 
large deviation functional in the form of using ([7]) and ([H]). (See the argument below 
(124]) in Sec. K2\ ) 

The physical argument developed above cannot be applied to non-equilibrium cases 
(po < Pi) because thermodynamics for non-equilibrium steady states has not yet been 
established. Nevertheless, it is expected that the large deviation functional for non- 
equilibrium cases might have a correspondence with the local equilibrium form, where 
the function F() given in is replaced with 

^ie(p()) = / dx[f(p(x)) - p(x)p(x)]. (9) 

The chemical potential in this expression depends on x, and its functional form is given 
by the local equilibrium thermodynamics. That is, 

p(x) = f'(p(x)). (10) 

Further, considering the relation 

D(p)d x p = a{p)d x f(p) (11) 
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derived from ([HI), we obtain 

J = -a{p)d x p. (12) 

This implies that a is indeed the conductivity. 

Here, it should be noted that the large deviation function is not equal to 
/3[Fi e (p()) — Fi e (p())], but there is a deviation from the local equilibrium part. We 
then express the function F in ([6]) as 

f(pO) = f 1c (pO) + n(pO)- (is) 

At present, the specific objective of this study is to find a physical interpretation of 
N(pQ) by obtaining its expression. 



3. Formal analysis 



3.1. Fundamental equality 
We first define 



Fo(p()) = / dxf(p(x)), 
Jo 

where / is given in dSj). We then rewrite as the continuity equation 
dtp + d x j = 0, 



(14) 



(15) 



with 



J 



-a{p)d x 



SF 

Sp(x] 



where we have used ( ITT]) . Since £ obeys the Gaussian process 
distribution of p( , t), < t < r, with fixed p( ,0), is written as 

SF 



V(pQ) ~ exp 



I 
4e 



/ dt I dx 



1 



j + a(p)d x - 



(16) 

the probability 



;i7) 



a(p) V ' Sp(x) 

where the current j is connected with p by the continuity equation given in f|T5l) . Note 
that the so-called Jacobian term does not appear in the leading order in the limit e — > 0. 

We next assume that the configuration p at £ = obeys the local equilibrium 
distribution 

P 



Ple(pQ) ~ eX P 



Let A(joQ) be an arbitrary function of pQ, and let (A) represent the average of A with 
respect to the path probability distribution P\ e (p)V(p()). 

In the equilibrium case (p — pi), the model satisfies the detailed balance condition. 
This property is also known as stochastic reversibility, and it can be formalized by 
introducing the time-reversed trajectory of p that is denoted by p' . This can be explicitly 
written as p'(x, t) = p(x, r — t). In terms of p* , the stochastic reversibility is expressed as 
a symmetry property (A) = (A^), where A^(pQ) = A(p^Q). The equality (A) = o4.t) is 
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not valid for non-equilibrium systems. We then attempt to derive its extended relation. 
The key equality for deriving this relation is as follows: 



Vip{)) - <'x P (--;/•„(/,( .;!) /•;,(/;( •())) 



(19) 



WO) *A e 

+ [ dt(p(l)j(l,t)-p(0)j(0,t)) 
Jo 

which is obtained by direct calculation. Using this, we calculate 
(A) = J VpPM ,Q))T{p{))A(p()) 

V P PMMV( P ()) ff\^f^ W ) 
P\ c {p{ ,0)) V(pQ) 

A f exp 0?- dt jf dx(d x p)j^ ^ . (20) 

Since the relation (A) = (A) is derived in the equilibrium case (c^/i = 0), (1201) 
is regarded as an extended form of the stochastic reversibility. It has been known 
that this equality yields many non-trivial relations, including the Green-Kubo relation, 
Kawasaki's non-linear response relation, and the fluctuation theorem [121 . 



3.2. Stationary distribution 

Following the method presented in Ref. [12] . we derive an expression for the stationary 
distribution P s (pQ) based on the relation in (1201) . The probability distribution of the 
configuration p is given by 

P s (p())= lim (6(p( ,r) - p())) ■ (21) 
Substituting A(pQ) = 8(p( , r) - p()) in ([20]), we obtain 

Ps(pO) = PM)) [ Vpe-&W>\ (22) 

J P ( ,o)=p() 



with 

roo pi r -r 

(23) 



/*00 /» 1 

S(pQ) = / dt dx 
Jo Jo 



1 -^^^(p)^) 2 -^^) 



Substituting (Q, fTTSl . and f|T8|) into fl22|) . and applying the limit e — > 0, we obtain 

N(P()) = J n min S(p()). (24) 

4 p();p( ,0)=p() 

In equilibrium cases, using the condition d x p = 0, we can trivially solve the minimization 
problem as N(pQ) = because S(p()) > for any p() and the solution of the diffusion 
equation with p( , 0) = p() yields E = 0. On the contrary, the minimization problem 
is not easily solved in non-equilibrium shown below. Note that the expression 

of the large deviation functional in terms of the minimization over trajectories was 
presented in Ref. |T3l [151 US] . However, to the best of our knowledge, the expression of 
N(pQ) given in (|2"4"]) has never been reported. 



Large deviation junctionals in fluctuating hydrodynamics 



7 



In order to determine the trajectory minimizing E on the right-hand side of 
we consider the variation p(x, t) — > p(x, t) + 5p(x, t). We first assume that the trajectory 
pQ that minimizes S(p()) satisfies 

lim p(x,t) = p(x). (25) 

t— too 

Further, defining a new variable 

d x u = (j + D(p)d x p)/a(p), (26) 
we calculate = + — as 

5S= / rft/ dx{[- ( r'(p)(9 x . M ) 2 -2 J D(p)(^ M )](5p) + 2(9 x . M )(5j)}. (27) 
Jo Jo 

From ffl5l) . 5j is related to <5p in the form 

6j(x,t) = -d t d x [ dyA- 1 (x,y)6p(y,t), (28) 



where A 1 (x, y) is the Green function satisfying c^A 1 (x, y) = 8(x—y), and A 1 (0, y) = 
A _1 (l,y) = 0. The substitution of (I28I) into the variational equation <5E = yields 

2 / dy(d y A-\y,x))d t d y u(y,t) = a\p){d x u) 2 + 2D{p)d 2 x u. (29) 
Jo 

To this point, u(0,t) and u(l,t) can be any function of time. We here assume that the 
trajectory minimizing E is obtained in a class of trajectories that satisfy 

u(0, t) = Uq and u(l,t) — Ui, (30) 
where Uq and U\ are constants in time. Then, (1291) becomes 

9 t « = - ( 7 / (p)(^w) 2 /2 - D(p)d 2 x u. (31) 
Furthermore, ffl5l) and (126]) lead to 

d t p = d x (D(p)d x p - a(p)d x u). (32) 

In this manner, with the two assumptions ( |25i) and ( 1301 . we have found that the 
trajectory minimizing E satisfies a set of equations (|3~T!) and (1321 . Here, it should 
be noted that the difference between (TT71) and (|23ll is the term jd x p in ( !23l) and that the 
variational equation for ( |23l) does not depend on this term. Thus, ( |3T1) and ( |32l) also 
determine the most probable trajectory under given boundary conditions in the limit 
e — > 0. Indeed, ( 13TI) and ( |32j) are identical to the Hamiltonian equation providing the 
most probable trajectories [13] . 

Using this property, we consider the minimization problem. First, let p F be the 
solution of the deterministic equation d t p = d x (D(p)d x p) with p(x, 0) = p(x). Then, 
(p, u) = (p F , const.) is one solution of the variational equation given in ( |3TT) and 
( 1321) . This solution describes the most probable relaxation behavior, provided that 
p( , 0) = pQ, but it yields S — > oo unless <9 x /i = 0. (See (|23|) .) Therefore, this most 
probable trajectory is not relevant in the minimization problem. 
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In order to seek another solution, we further assume that the trajectory minimizing 
E satisfies 

lim u(x, t) — u*(x). (33) 

Then, j approaches in this limit, where 

j* = -D(p)d x p + a(p)d x u*. (34) 

Substituting this expression into ( 13T1) . we obtain = ±J. When = J, we derive 
d x u* = 0, which corresponds to the case p = p F . We thus pay attention to the solution 
for the case j* = — J. Then, ( 134"1) leads to 

d x u*(x) = 2d x p, (35) 

and we can easily confirm that S < oo for this solution. 

The physical interpretation of this solution is obtained by considering the equality 
j* = —J. Let p B (x,t) = p(x,—t) for the solution. Then, the current for p B is J at 
t = — oo. Further, defining 

u B (x, t) = -u(x, -t) + 2f'(p(x, -*)), (36) 

we can directly confirm that (p B ,w B ) satisfies fl3T|) and fl32l) with the condition 
(p B ,w B ) — >■ (p(x),0) as t — > — oo. Thus, p B (x, t) describes the most probable process of 
generating the fluctuation p() at t = starting from p(x) at t = — oo. Note that the 
boundary values w B (0,i) = u B (l,t) = are consistent with the condition imposed in 

Ref. rng . 



4. Perturbation theory 

4-1. Idea 

Before developing a perturbation theory, we consider a numerical method by which we 
can calculate N(pQ). First of all, we must solve fl3T|) and fl32|) under the conditions 
(p,u) — > (p, ii*) as £ — > oo and p( , 0) = p(). Note that u(x, 0) should be determined 
such that (p, m) — > (p, m*) as £ — > oo. However, it is quite difficult to perform a shooting 
method for partial differential equations. 

Instead of such a direct method, we devise the following recursive procedure. 
Initially, we set 

m (0) (x,£) =u*(x). (37) 
Then, for n > 0, we numerically solve 

dtpM = d x (D{pW)d xP W - <r(pV>)d x uW) (38) 
with 0) = p(x), and next we consider 

d t u {n+1) = -a'(p {n) )(d x u {n+1) ) 2 /2 - D(p {n) )d 2 x u {n+1) , (39) 
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with u( n+1 \x, oo) = u*(x). Here, solving (I3"9~j) . we set tt(x, — £) = u(x,t). Then, u 
satisfies 

d t u {n+1) = a f (p {n) )(d x u {n+1) ) 2 /2 + D(p^)<9;^ (n+1) , (40) 

for — oo < t < 0. Because tt^ n ^(x, — oo) = w*(x), we can have a standard algorithm of 
the forward time evolution starting from u*(x), and as the result vS n \x, 0) = u^ n \x, 0) 
is determined. If we confirm the convergence p( n '(x,t) — > p(x,t) and u^ n \x,t) — > u(x,t) 
within the desired numerical accuracy, we obtain an approximate solution for (I3T1) and 
( 1321) under the conditions (p, w) — > (p, u*) as t — > oo and p( , 0) = pQ. 



Perturbative expansion 

The method we considered above is also useful for developing a perturbation theory. 
As the simplest example, we derive N(pQ) up to the order of J 2 . Checking the J 
dependence of N given in ( 1241) . we find that the 0(J 2 ) contribution can be ignored in 
the estimation of u and j. Considering this, we investigate the recursion equation. First, 
we have u^°\x,t) = u*{x) and 

d t pW = d x (D(p^)d x p^ - a{p^)d x u*) (41) 

with p^(x, 0) = p(x). Next, we solve 

d tU W = -D(p^)d 2 x u^ - a\p^){d x u^f/2 (42) 

under the condition lim^oo u^(x, oo) = u*(x). Substituting u^(x, t) = u^°\x,t) + 
v^(x,t) into this equation and expanding in terms of J, we obtain 

v {1) (x,t) = 0(J 2 ). (43) 

Since u^(x,t) = + 0(J 2 ), we have arrived at the fixed point solution by ignoring 
the 0(J 2 ) contribution in u and p. Substituting this solution into (I23l) . we obtain the 
expression of iV(p()) as 

N(p())=(3 dt / dx(d x p)[D(p^)d x p^-a(p^)d xf i}+0(J 3 ), (44) 
Jo Jo 

with flUD and p^(x,0) = p{x). 

Now, let us recall the physical interpretation of the solution. As discussed in the 

last paragraph of Sec. [3j p B (x,t) = p^(x, —t) describes the most probable process for 

generating the fluctuation p() at t — starting from p(x) at t = — oo. Further, from 

fj4T|) . it is found that J B = D(p B )d x p B — a(p B )d x ii* represents the particle current in 

this most probable process. Using this current, we rewrite N(pQ) as 

N(f>0) = -P f dt [ dx(d xf i)[-J B - o-(p B )d xf i] + 0{ J 3 ). (45) 

J -oo Jo 

Here, —f3(d x p)J B is the entropy production rate observed in the process, while 
f3a(p B )(d x p) 2 is interpreted as the entropy production rate under the assumption 
that the steady state with the density p B (x,t) is locally realized in a space-time 
point (x,t). We call the latter quasi-steady entropy production. We then interpret 
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(3(d x fi)[— J B — cr(p B )d x p] as an excess entropy production ratio. With this interpretation, 
we can regard N(pQ) as the space-time integration of the excess entropy absorption rate 
during the most probable process of generating the fluctuation pQ at t = starting from 
p(x) at t = — oo. (Note the minus sign before the space-time integration.) 



4-3. Small fluctuations 



When we focus on small fluctuations, we can make another approximation by which a 
quantitative calculation becomes possible. Indeed, substituting p(x, t) = p{x) + <p(x, t) 
into dl!]) , we obtain 

J 2 A(p ) " < l 



N(p()) 



2a(p y 



dt / dx((f)(x,t)) 2 + O(J 3 ,0 3 ) 
n Jo 



(46) 



with 



A(p ) = a"( Po ) - a'(p )D'(p )/D(p ). (47) 

Furthermore, since we ignore the terms of 0(J 3 ) in N, we can assume that satisfies 
d t 4> = d x (D(po)d x 4>) with the initial condition 4>(x,0) = p(x,t) — p and the boundary 
conditions 0(0, t) = 0(1, t) = 0. Here, considering that p 1 — p ~ O(J), we find that p 
in ( 1461) and the boundary conditions can be replaced with an arbitrary value in [p ,Pi]- 
This approximation leads to 

(4f 



P B (pQ) ^ exp 
Here, we have defined 

D(p(y)) 



L(y,y') 



§■ f dy I dy'L{y,y')4>{y)4>{y' 
Ae Jo Jo 

J 2 A(p ) 



Kv-y') 



&(p(y)) °(po) 2 jo jo 



oo rl 

dt I dxG(x,y,t)G(x,y',t),(<i9) 



where G(x, y, t) is the Green function that satisfies 

[d t -D(p o )dl]G(x,y,t) = (50) 

when t > 0, G(x, y, 0) = 5(x — y), and G(0, y, t) = G(l, y, t) = 0. An expression that is 
essentially the same as that in (T49l was derived by the analysis of the linearized equation 
around the steady solution p [7]. 

Furthermore, note that the space-time integration appearing in fT49|) is calculated 



as 

dt [ dxG(x, y, t)G(x, y', t) = V -^-Ae{y - y') + (y <- y'), (51) 
o Jo ^ u \Po) 

where 8( ) is Heaviside's step function. Then, from the Gaussian nature of the 

fluctuations, we derive 

(0(2/ W)> 



+ eT 



D(p(y)) 
J 2 A(p ) 



D(pof 



y'(i-y) 



o(y - y') + (y^ y') 



(52) 
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In order to present a simpler demonstration, we calculate the intensity of the fluctuations 
of the spatially averaged density 



The first term on the right-hand side represents the contribution of local equilibrium 
fluctuations, and the second term originates from non-local fluctuations. For the special 
case D = 1, a = p(l — p), and T = 1, the second term becomes — (pi — po) 2 /12, which 
is consistent with the result obtained from the exact solution in Ref. [I]. 

5. Concluding remarks 

We have derived the large deviation functional for a simple model using fluctuating 
hydrodynamics. In particular, focusing on the deviation N(pQ) from the local 
equilibrium part, we obtain its lowest order expression with respect to the average 
current J. This expression provides us with the physical interpretation that N(pQ) 
corresponds to the excess entropy absorption during the most probable process of 
generating the fluctuation p() at t = starting from p(x) at t = — oo. 

The evolution equation describing the most probable process generating the 
fluctuation p() at t = is called the adjoint hydrodynamic equation. Within our 
approximation, we have obtained the adjoint hydrodynamic equation as ( 14T1) . In contrast 
to equilibrium cases, the generating process is not given by the time reversal of a 
relaxation process owing to the existence of the second term in (I4T1) . Rather, it may 
be noticed that (jUJ) is similar to the equation for driven diffusion systems. It is an 
interesting subject to elucidate the physical picture for this asymmetry. 

Our results can be generalized to those for several cases in a straightforward 
manner. For example, the calculation of the large deviation functionals for driven 
diffusive systems may be a natural problem, which will be studied in a similar manner. 
Furthermore, the analysis of fluctuating hydrodynamics for a simple fluid is a physically 
important problem [T5], [T7] . With regard to this problem, it has recently been found that 
the stationary distribution for non-equilibrium systems with multiple heat reservoirs is 
expressed in terms of the excess heat up to the order of the square of the average heat 
flux [18]. The connection between the two results will be explored. 

Before concluding this paper, let us recall again that large deviation functionals in 
equilibrium cases are expressed by a thermodynamic function. Thus, we are naturally 
led to consider a thermodynamic framework consistent with the expression of the large 
deviation functional that we have calculated in this paper. For example, by using 
a method similar to that developed in Ref. [TH], we can derive an identity for the 
operations expressed by a parameter change. Applying the identity to the case in which 




(53) 



Using (1521) . we obtain 




(54) 
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the boundary conditions are changed, one may find some insight related to an extended 
framework of thermodynamics. 

Furthermore, we wish to highlight the fact that the additivity of the system is 
an important property in a thermodynamic framework [3j. However, the additivity 
at a thermodynamic level is not directly related to the additivity of large deviation 
functionals in non-equilibrium steady states because of the existence of non-local 
fluctuations. Therefore, one may study the additivity principle that was discovered 
for exactly solvable models [U EJ E] from the viewpoint of an extended framework of 
thermodynamics. The concrete expression given in (j4"5]) will help us to consider such a 
problem. 
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